Novel insights into the 77* — > n° transition form factor 
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BaBar's observation of significant deviations of the pion transition form factor (TFF) from the 
asymptotic expectation with Q 2 > 9 GeV 2 has brought a serious crisis to a fundamental picture 
established for such a simplest qq system by perturbative QCD, i.e. the dominance of collinear 
factorization at high momentum transfers for the pion TFF. We show that non-factorizable contri- 
butions due to open flavors in 77* — > ty° could be an important source that contaminates the pQCD 
asymptotic limit and causes such deviations with Q 2 > 9 GeV 2 . Within an effective Lagrangian 
approach, the non-factorizable amplitudes can be related to intermediate hadron loops, i.e. K^*' 
and etc, and their corrections to the ir° and 77 TFFs can be estimated. 

PACS numbers: 13.40.Gp, 12.38.Lg, 11.40.Ha 



I. INTRODUCTION 

Since the foundation in late 1970s perturbative 
QCD (pQCD) has been a powerful tool to explore the 
strong interaction phenomena at large momentum trans- 
fers. This is such a kinematic region where the exclusive 
transition matrix element can be expressed as the convo- 
lution of the perturbative calculable coefficient and non- 
perturbative light-cone distribution amplitude (DA) that 
describes the longitudinal momentum fraction of quarks, 
namely the so-called collinear factorization. In particu- 
lar, with the availability of clean electromagnetic probes, 
the photon-pion transition form factor (TFF) will be the 
most ideal subject to test the validity of pQCD and par- 
tonic structure of pion. 

Experimentally, the pseudoscalar mesons' TFFs have 
been measured by CELLO and CLEO collaborations 
with the virtuality of one photon Q 2 up to 9 GeV 2 and 
the other nearly on-shell @, Q ■ The data fit very well 
the pQCD leading order interpolation formula Q , which 
stitches together the chiral anomaly and pQCD asymp- 
totic limit, 
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where f„ = 92.3 MeV, f v = 97.5 MeV are the cor- 
responding decay constants of the pseudoscalar mesons. 
In fact, it has been shown that nearly all the analyses up 
to next-to-leading order (NLO) radiative and twist-four 
corrections favor the endpoint-supprcsscd or asymptotic- 
like pion DA, which means that the first two nontrivial 
Gegenbauer moments of the DA are small and have op- 
posite signs in the region of Q 2 = 9 GeV 2 [H, @ . 

However, the celebrated confirmation of pQCD dom- 
inance at large momentum transfers was shortly over- 
whelmed by the surprising BaBar measurement 0, H| . It 



shows that the rescaled pion TFF exhibits a continuous 
growth at 9 < Q 2 < 40 GeV 2 , even beyond the asymp- 
totic limit of pQCD. In contrast, the 77 and rj' TFFs seem 
to be consistent with the pQCD expectations. Thus, a 
coherent understanding of the tt° and 77 (77') TFFs would 
be a necessity for resolving this puzzle. 

The BaBar observation immediately arouses tremen- 
dous interests in the pion TFF. In the pQCD scheme, 
the scaling violation of the BaBar data can be better ex- 
plained by either a wider DA with deep midpoint or even 
a flat one [9T4l3j. However, they both conflict with some 
other constraints such as the data at low Q 2 , midpoint 
value of the pion DA from light-cone sum rules [14j . or 
operator product expansion_(OPE) analysis for the lead- 
ing handbag contribution [9| etc. Moreover, a flat DA 
means pion is an elementary field and can interact with 
quarks locally. Thus, a quark loop can also give the log- 
arithmic enhancement but with rather low constituent 
quark masses [l~5T - tl7| . In this case, the QCD evolution 
is switched off, and the whole process is totally non- 
factorizable 911. Therefore, a rather flat DA with van- 
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ishing endpoints may still be a temporary solution 14 1 
taking into account that it is not supported by the u, d 
quark behavior in the r\ and 77' TFFs [3. It is even 
claimed that a peculiar mechanism beyond the standard 
model is needed to explain the BaBar puzzle fl9l - l2lj |. 

In the chiral limit approximation, there is only one 
large physical scale in the hard part of pQCD calcula- 
tion. So the power law is explicit and concrete. In turn, 
the abnormal growth of the rescaled pion TFF implies 
that an additional mass scale would be present. A nat- 
ural mechanism for generating a mass scale should sat- 
isfy the following constraints: i) It should keep the gen- 
eral OPE analysis. In another word, the solution should 
originate from non-OPE or non-power corrections if only 
light quarks are involved; ii) The new mechanism should 
leave the chiral anomaly intact at Q 2 = since the room 
for the theoretical and experimental improvement is very 
small [22I]; iii) The pQCD asymptotic prediction should 
be valid, though the current experimental region in the 
exclusive processes seems to be still far away from the 
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asymptotic region determined in inclusive processes [23| ; 
iv) The flavor symmetry breaking of the DA should be 
explained in the same scheme, and it is better to keep 
SU(3) flavor symmetry taking into account its big suc- 
cess in the effective theory; v) Pion should be treated 
as a bound state of current quarks. It is necessary to 
make the new mechanism compatible with the numerous 
analyses for both the asymptotic-like DA and collinear 
factorization before the BaBar data since they have been 
cross checked based on different powerful tools. 

To find a solution that satisfies the above constraints 
simultaneously seems not easy. In the literature there are 
also proposals to suggest that the BaBar puzzle may need 
explanations beyond QCD. However, before jump to any 
unusual solution, we should first have a full understand- 
ing of the evolution of non-perturbative mechanism from 
low to high Q 2 . In principle, any solution satisfying those 
constraints based on QCD should be treated seriously. 

In this work, we propose a novel prescription to address 
above questions as a compensation of pQCD. Since non- 
perturbative flavor changing processes (FCP) may play 
a role as "mass transmutation" above certain mo- 
mentum transfers, we assume that such non-perturbative 
mechanisms can be expressed in terms of hadronic de- 
grees of freedom. With an effective Lagrangian approach 
(EL A) , we shall show that intermediate meson loops have 
a peculiar evolution in terms of Q 2 , which could be an 
important mechanism causing the TFF deviations from 
the asymptotic behavior in the region of 9 < Q 2 < 40 
GcV 2 . 

BaBar Collaboration recently also measured the r] c 
TFF (24|. A continuous growth of the rescaled TFF 
(or slow decrease of the TFF) is observed clearly with 
the space-like momentum transfer square up to 40 GeV , 
which is caused by the large charm quark mass (25j . It 
is also shown that a small charm component in 7/ , which 
is needed to explain the abnormally large branch ratio of 
B — > Kr( and B — > X s r]' [2(| l27j . may cause large de- 
viations from the asymptotic prediction in the medium 
momentum transfer region while affects little the low mo- 
mentum transfer region [28| . So it is natural to speculate 
that the charm quark may play a role in the deviations 
from the asymptotic predictions in the BaBar kinematic 
region. The first possible QCD diagram is the photons 
fuse into a charm quark loop which then couples to two 
gluons. These two gluons then evolve to the pion through 
triangle anomaly. This factorizable charm loop is also a 
possible correction for the chiral anomaly (22|, and can- 
not be recognized by the OPE based only on the massless 
quarks. Although the correction at Q 2 = 0, i.e. the chi- 
ral anomaly, is tiny, its evolution with the photon virtu- 
ally could make it relatively important in larger photon 
virtuality. Unfortunately, for the case of pion transition 
it is an isospin suppressed process proportional to the 
u, d quark mass difference. Therefore, its influence on 77 
is much stronger than that on the pion. To avoid the 
isospin suppression, the u, d quarks would interact with 
photons more directly, and the gluon exchanges between 



charm quark and light quark must be non-perturbative 
to evade the strong coupling suppressions. Such a sce- 
nario can be equivalent to the introduction of the inter- 
mediate meson loops based on the quark-hadron duality 
argument. It will be our focus in this work to investigate 
such a non-perturbative effect on the pion TFF in the 
BaBar kinematic region. 

The rest of this paper is organized as follows: The ef- 
fective Lagrangian approach for the intermediate meson 
loops is introduced in Sec. II. Numerical results and dis- 
cussions are presented in Sec. Ill, and a brief summary 
is given in Sec. IV. 



II. EFFECTIVE LAGRANGIAN APPROACH 
FOR THE INTERMEDIATE MESON LOOPS 

The effective Lagrangian for the mesons 
(£)(*)0 ) _D(*)+ ) £)W+) couplings to light pseudoscalar 
mesons (2{| has the following expression: 
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where d a =d a — d a , and V denotes the pseudoscalar 
octet mesons 
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The corresponding Lagrangians for the photon and £)W 
couplings are 
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where = d^A, - d v A^ £>*«>+ = - d„Df>+, 

and the vertices D°D°j and D*°D*°7 do not exist. 

The coupling constants 3d»+_d+ 7 and gD*°D°i ar e ex- 
tracted from [3(j| : 



g D , +D+1 = -0.5 GeV" 1 , g D .o DO „ = 2.0 GeV" 1 , 
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where the signs and relative size can also be deduced 
from constituent quark model, see also [3l| for the lattice 
calculation. Similarly, the coupling constant gwov = 
QD'Dtt = 17.9 [IH, and gD'D*v can be obtained from 
the heavy quark effective theory 
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Four types of loops would contribute to the pion TFF, 
i.e. D*D(D), D*D(D*), D*D*(D*), and D*D*(D), 
where the mesons in the parentheses are the ex- 

changed mesons between two photons. The kinematic 
conventions arc illustrated by Fig. [TJ The transition am- 
plitude can be expressed as follows: 
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where I3 4^ are the vertex functions given by the effec- 
tive Lagrangians, and 0^=3,4.5 are the denominators of 
propagators of the intermediate mesons, respectively. 

At hadronic level, the calculation of intermediate me- 
son loops only includes S- wave ground state mesons. In 
our case, the intermediate charmed mesons include D 
(D s ) and D* (D*), and in the strange sector K and K* 
are considered. One essential question which is very of- 
ten raised is the contributions from higher excited meson 
loops apart from those S-wave states. Intuitively, one 
would expect that a complete set of intermediate meson 
loops should be included based on the quark-hadron du- 
ality argument. In our approach, the following reasons 
would allow us to only consider contributions from the 
S'-wave states. Firstly, for a given kinematic condition, 
we would expect that the local EM dipole couplings in 
Eq. ([7]) would be much larger than the excited states 
due to quantum mechanical selection rules. For electric 
dipole couplings, e.g. jD + D~. the ground state would 
be less suppressed by the wavefunction overlaps than ex- 
cited states. Secondly, in many cases the quark-hadron 
duality is locally broken such that a subset of hadron 
loops would play a dominant role in observables. More- 
over, we find in numerical simulations that apart from 
the vertex couplings the loop integrals general decrease 
when the masses of the intermediate mesons increase. In 
this sense the contributions from the ground state S'-wave 
intermediate meson loops can be regarded as a reason- 
able approximation. Meanwhile, although we still lack a 
systematic evaluation of how good such an approxima- 
tion would be in different circumstances, it interests us 
to explore such a possible mechanism that accounts for 
the non-factorizable contributions to the pion TFF. 

We also note that in the process of 77* — > 7r°, there ex- 
ists a significant difference between go*+D+ 7 an d to* B°7 
as given in Eq. ([7]). In some other diagrams, the neutral 
meson loops even vanishes, e.g. due to the 'yD°D° ver- 
tex. It means that the isospin violation effects are actu- 
ally enhanced by the electromagnetic interaction. Also, 
a nonvanishing charmed meson loop contribution must 
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FIG. 1: D'*' meson loop contributions to pion TFF. The 
other diagrams can be obtained by exchanging the photon 
momenta. 



be present in the low Q 2 region which comes from the 
destructive sum between the charged and neutral 
meson loops. At sufficiently high Q 2 the meson loop 
contributions should become negligible as expected by 
pQCD. The interesting question is at which range of Q 2 
such an effect becomes negligible. Meanwhile, an asso- 
ciated difficulty is to model-independently quantify the 
loop contributions in a broad range of Q > 0. In this 
sense our study is motivated to provide a qualitative es- 
timate of the intermediate meson loop effects as part of 
the non-factorizable contributions in terms of Q 2 . 

Generally speaking, the loops suffer from ultraviolet 
divergence which means that we should replace the lo- 
cal coupling constants by the non-local form factors to 
suppress the divergence. As usually done, we introduce 
the form factor F(p 2 ) to take into account the off-shell ef- 
fect of each vertex, which also reflects the substructure of 
external hadrons. With the large space-like momentum 
transfers, each vertex would be far off-shell. Therefore, 
we adopt the following tri-monopolc form in the calcula- 
tion: 
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This turns out to be a reasonable consideration by our 
numerical simulation. The cutoff A^ can be parameter- 
ized as A, = rrii + uAqcd with Aqcd — 220 MeV, and 
nii is the exchanged meson mass. Parameter a is usually 
taken at 0(1) and not necessarily the same for different 
open flavors. 

Apart from the D loops, the Kaon and B loops 
can also contribute. The Kaon loops are similar to 
Fig. [T] where the coupling constants can be deter- 
mined by experimental data [33|, i.e. gx'+K+-y = 



0.84 GeV" 1 , g K ->o K o~ 



-1.27 GcV-\ g K . +K+7r o 



4 



— gK*°K°Tv° = — 9.17/V2) .90i)7 = 0.69 GeV -1 and 
<707j'7 = —0.73 GeV -1 . The SU(3)f symmetry then 
gives gK* a K*°TT° = -9k*+k*+tt° = — ffp°um°/2 = 
-11.58/2 GeV -1 . 



III. RESULTS AND DISCUSSIONS 

The introduction of an empirical form factor will bring 
in model-dependent features. In order to give a meaning- 
ful interpretation of the results, we outline the conditions 
for constraining the parameter a as follows: 

i) For an adopted value of a, the meson loop correc- 
tions to the TFF at Q 2 = should be negligibly small 
in order not to conflict with the well established theoret- 
ical and experimental results for the chiral anomaly [22j . 
In practise, the D or Kaon meson loop contributions to 
the anomaly are restricted to be less than 0.005 GeV -1 
(absolute value) in comparison with Eq. ([T]), F 77T (0) = 
0.27 GeV -1 . This serves as a strong constraint for the 
upper limit of a. 

ii) The form factor of Eq. (|10[) introduces additional 



singularities into the integrals empirically [34| . which 
means that in order to reduce the model-dependence, 
parameter a should have a sufficiently large value, e.g. 
a > 1. This condition can be satisfied in 77* — > tt° since 
all the internal exchanged mesons are highly off-shell. 

iii) We neglect B meson loops because the numerical 
calculation shows that apart from couplings, their con- 
tributions are usually one order magnitude smaller than 
D and K loops. In addition, we do not include the non- 
strange light meson loops to avoid double counting. 

In fact, we find that the above conditions can indeed 
provide a stringent constraint on a, and 1 < oik < 2, 
1 < old < 4 can be determined in Fig. [31(b). 

In Fig. [2] we plot pion TFF given by the exclusive D 
and Kaon loops with ud = 3 and uk — 1.5 in terms 
of Q 2 . It should be stressed that, the meson loop con- 
tributions to the chiral anomaly are nearly zero because 
of the destructive interferences from the D and K loops. 
It is indicated by the thick solid line at zero momentum 
transfer point in Fig. [2] The dominance of D loops over 
the Kaon loops can be recognized. In particular, the 
D*D*(D) loop exhibits an enhancement in the BaBar 
kinematics due to the relatively large coupling differences 
between the charged and neutral meson loops as shown 
in Eq. ([7]), i.e. the vertex coupling product of the neu- 
tral D meson loops is 16 times larger than that of the 
charged D mesons. The Kaon loops are much smaller 
because of the destructive interferences among different 
loop amplitudes and smaller oik- 

A sensible feature of the D loops is that, the total me- 
son loop contribution decreases faster than pQCD with 
the increasing Q 2 , Thus, it guarantees the asymptotic 
prediction of pQCD at high Q 2 . This is different from 
the logarithmic [T(| or double logarithmic [HI Q3] growth 
predicted by some analyses with fiat DA or equivalent 
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FIG. 2: (a) D meson loop contributions to the pion TFF. 
The thin solid, dashed, dot-dashed, and dotted lines denote 
D*D(D), D*D(D*), D*D*(D*), D* D* (D) loops, and the 
thick solid line for the total loops, dot-dot-dashed for Eq. {T}. 
(b) Kaon loop contributions with the same notations as (a), 
but replaced by Kaons. 
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FIG. 3: (a) The rescaled pion TFF in comparison with the 
experimental data. The dashed line is given by the pQCD 
interpolation formula. The meson loops combining the pQCD 
results give the shadowed band with cxd = 2 and an — 3 for 
the lower and upper bound, (b) The a dependence of the 
meson loop contributions to the chiral anomaly. The dashed, 
dotted, solid lines are for D loops, Kaon loops, and their sum 
with the same a, respectively. 



pion vertex function, although we should caution that the 
EL A result at high Q 2 is inevitably model-dependent. 

In Fig. [3Ja), the shadowed band is the exclusive con- 
tributions from the intermediate meson loops plus the 
pQCD prediction, which in this paper means the interpo- 
lation result of Eq. ([T]). It should be noted that there is no 
double counting between our meson loops and the pQCD 
interpolation result involving only the u, d quarks. The 
range of the shadowed area is given by an = 2 (lower) 
and 3 (upper bound) with ax = 1-5 fixed. This is a 
rather conservative estimate of the meson loop contribu- 
tions based on the above constraints on a. It shows that 
the inclusion of the meson loops significantly improves 
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FIG. 4: (a) The rescaled 77 TFF in comparison with the exper- 
imental data. The thick line denotes meson loops plus pQCD 
interpolation result which is the dashed line, while the dot- 
dashed, dotted and thin solid lines are for the \nn), \ss), and 
physical 77 from meson loops, (b) The a dependence of the 
meson loop contributions to the chiral anomaly. The dashed, 
dotted and solid lines are from the D loops, Kaon plus (ficfiiv 1 ''' 1 ) 
loops, and their sum with the same a, respectively. 



the region of Q 2 > 9 GeV 2 while affects little in Q 2 < 5 
GeV 2 . Therefore, the CELLO, CLEO and BaBar exper- 
iments can be explained simultaneously without modify- 
ing conclusions obtained before. We also note that the 
discrepancies with the data at low Q 2 are mainly caused 
by the pQCD interpolation formula Eq. (TTJ) which would 
break down due to the growing importance of the QCD 
corrections with negative sign @. We also emphasize 
that we do not try to fit the data and judge the possible 
configuration of the pion DA. The meson loop mechanism 
only tells us that in the high momentum transfer region, 
the asymptotic-like pion DA may be still compatible with 
the BaBar data. 

The meson loops can also contribute to the 77 and 77' 
TFFs. It is thus natural to expect that the same formula- 
tion in the SU(3) f symmetry limit should keep consistent 
with the theoretical and experimental results for 77 and 
77', or at least for ?]. Interestingly, the data @ show that 
the rise of the rescaled TFFs is much weaker than the 
rise of ir° case. It implies that the u, d components in 77 
and 77' would be significantly different from that in ir° if 
the analysis is only based on pQCD Therefore, it is 
interesting to examine the meson loop contributions here 
as a direct test of the proposed mechanism. It should be 
pointed out that a strict anomaly sum rule shows that 
the same non-perturbative non-OPE correction to the 
continuum as in t he p ion TFF should also be present in 
the 77 and 77' TFF (^H, although the 77 TFF favors the 
result from the pQCD interpolation formula. Here we 
do not discuss the 77' TFF because of its possible mixing 
with the gluonic and/or charm component (37j . 

In Fig. 21 the rescaled 77 TFF is plotted. Considering 
the 77-77' mixing, an extended interpolation formula for 77 



and 77' TFF was given by Foldmann and Kroll [3£ 
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where C q = 5/9^2 and C s = f/9, and f q and f s arc 
decay constants of the qq = (qq + dd) / \[2 and ss compo- 
nents. The following couplings are extracted in Ref. [39|: 
f q = 1.07/,,, f s = 1.34/,, with f n = 0.131 GeV, and the 
coupling transformation satisfies: 



fq fs 
J rj J rj 

j n' J rj 1 



U(a P ) 



f g 
/, 



(12) 



where the nonstrange qq and ss mixing angle favors a 
range of a P = 6 P + 54.7° ~ 39° - 43° @, M, M, EH- 
With the above coupling constants and flavor mixing an- 
gle ap = 39.3° given by Ref. [3!|, one can see clearly that 
the result of Eq. (TPTj) does not distinguish significantly 
from that of Eq. ([1]). Remember that our intermediate 
meson loop contributions do not interfere with the pQCD 
interpolation result. A further delicate treatment of the 
interpolation formula can be regarded as an improvement 
from the perturbative side, while our intermediate meson 
loops are from the non-perturbative side. 

In the meson loops the corresponding couplings (in- 
cluding D s loop) are gD*Dqq(q = u,d,s) = go*DV and 
go* D»qq extracted from Eq. ©. With SU(3) f symmetry, 
there are g K *K q q = -9.17, g K -K*qq = 11.58/\/2 GeV" 1 , 
and g^ss — 11.58 x y/2 GeV -1 for the additional (t><j>{r]^ ) 
loops. Including the 77-77' mixing angle, the 77 couplings 
to different mesons can be obtained: go*+D+ v = 

9D*Dqq cos ap/v2, gD'°D° v = gD'DqqCosa P /\/2, and 
gD*D s jj = -gD'DqqSinap. 

The same form factor parameters, ap> = 2.5 and ax = 
cf^ = 1.5, as in the range of 7T° case are adopted. Note 
that in Fig. HJb), the constraints of a are also satisfied 
by our setting. The total meson loop correction is much 
weaker than in the pion case since the qq = (qq + dd) / \/2 
and ss components significantly cancel out in Q 2 < 40 
GeV 2 . Our result suggests that the BaBar data can be 
well understood without jeopardizing those old analyses 
of the 77 TFF. 

The datum at Q 2 = 112 GeV 2 is actually measured at 
the time- like point Q 2 = -112 GeV 2 [if. Due to the 
analyticity of QCD, the time and space-like form factors 
can be related to each other at Q 2 — > 00. We mark the 
time-like datum point in Fig. [4] assuming that the \Q 2 \ 
value for the time-like virtual photon is high enough for 
the analytical continuation of the TFF. In this sense, it 
may provide a guidance for the magnitude of the corre- 
sponding space-like form factor at Q 2 ~ 112 GeV 2 . It is 
interesting to see that our result converges at high Q 2 as 
expected. Namely, the meson loop contributions should 
vanish at high Q 2 . Although the pQCD interpolation 
formulas of 77 and 77' are not as well-established as that of 
pion, e.g. radiative and higher-twist corrections need to 
be systematically included, the role played by the meson 
loops turns out to fit well the observed pattern. 



6 



IV. SUMMARY 

In summary, we have proposed that the meson loops as 
a non-pcrturbative component of the pion TFF may still 
play an important role up to Q 2 ~ 40 GeV 2 , hence cause 
deviations from the collinear factorization results. This 
mechanism is non-factorizable and different from the fac- 
torizable charm quark loop at quark-gluon level. The 
factorizable charm quark loop in the pion TFF must be 
suppressed by the u-d quark mass difference as compared 
with the 77. The meson loops in terms of quark-hadron 
duality may correspond to some continuum corrections in 
the spirits of the strict anomaly sum rule (35j . It should 
be stressed that this solution keeps consistent with the 



chiral anomaly without violating the SU(3) f symmetry. 
Nevertheless, this mechanism does not bring conflicts to 
against the pQCD analyses before the BaBar result, and 
is an economic explanation for the BaBar puzzle. 
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